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dimensional  flows.  The  other  involved  a  streamline-diffusion  finite  element 
method  for  inviscid  flows. 
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FINAL  REPORT 

"Computational  and  Experimental  Fluid  Dynamics" 

(MIR  Grant' for  the  period  June  1984  through  September  1985 

m * 

Funds  were  used  for  two  major  categories:  support  for 
released  time  for  the  principal  investigator  and  salary  for  a 
programmer.  The  former  allowed  the  principal  investigator  to 
visit  both  Los  Alamos  National  Laboratory  and  the  Mathematics 
Research  Center  at  the  University  of  Wisconsin  during  the  winter 
and  spring  of  1985,  as  well  as  providing  50%  released  time  during 
the  fall  of  1984.  This  facilitated  the  research,  especially  with 
regard  to  collaboration  with  W.  G.  Pritchard,  who  was  visiting 
the  Mathematics  Research  Center  during  that  time .  Research  on 
both  internal  and  free-surf ace  flows  related  to  experimental 
situations  was  continued.  Part  of  this  research  was  based  on 
using  a  low-order  code  developed  previously,  jointly  by  Pritchard 
and  Scott.  Funds  for  the  programmer  allowed  two  substantial  new 
programming  projects  to  be  initiated  during  this  time.  One 
concerned  investigation  of  techniques  for  implementing  a  quartic 
finite  element  method  for  viaooua,  incompressible,  two-dimensional 
flows.  The  other  involved  a  streamline-diffusion  finite  element 
method  for  inviscid  flows. 

The  collaboration  with  Pritchard  during  this  time  period 
focused  on  two  experimental  situations.  One  involves  internal 
flow  in  a  complicated  geostt ry  (see  Figure  1,  which  includes  an 
indication  of  a  typical  mesh  used  in  computations)  that  models  the 
critical  part  of  a  so-called  "stress  meter"  for  determining 
rheological  properties  of  non-Newtonian  fluids.  Experimental  work 
on  this  was  done  by  Pritchard  and  Arthur  Lodge  of  the  University 
of  Wisconsin.  A  pleasant  result  of  this  part  of  the  research 
during  this  period  was  that  we  discovered  a  nanuf aoturing  flaw  in 
a  new  meter  by  using  the  numerical  code  to  predict  a  particular 
stress  difference.  A  second  result  of  this  research  is  a  series 
of  comparisons  of  different  numerical  methods,  done  jointly  by 
David  Malkus  and  John  8trikwerda  of  the  University  of  Wisconsin 
together  with  Pritchard  and  Scott.  It  is  hoped  that  this  will 
help  to  establish  the  "stress  meter"  model  problem  as  a  benchmark 
with  a  firmer  physical  foundation  than  the  f requently-used  "driven 
cavity"  problem. 

The  typical  quantity  of  interest  in  the  problem  depicted  in 
Figure  1  is  the  normal  streas  difference  between  the  bottom  of  the 
slot,  or  hole,  in  the  middle  of  the  figure  and  the  top  of  the 
channel  above  the  hole.  This  differnce  is  sero  for  Reynolds 
number  equal  to  sero,  and  it  is  plotted  as  a  function  of  Reynolds 
number  in  Figures  la  and  lb  for  different  hole  depths  (Figure  1 
depicts  a  hole  depth  of  1.0).  Note  that  the  stress  difference 

wnnr'tnni«*al  I  W  ,  aa  a  f  nnr»  i  r»n  for  dept.h«  Hoim 

to  0.75.  however,  for  *  depth  of  0.5  it  suddenly  jumps  up.  This 
needs  to  be  studied  in  aore  detail. 
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The  other  experimental  phenomena  under  study  during  thia 
period  involved  free -surf ace  flow.  Pritchard  had  recently 
completed  a  series  of  experiments  at  the  University  of  Essex,  and 
work  was  begun  to  modify  the  code  in  order  to  allow  for  a  free 
surface.  Figure  2  shows  a  typical  plot  of  the  experimental  data; 
note  the  complexity  of  the  free  surface,  even  though  this  is  a 
relatively  low  Reynolds  number.  Our  goal  will  then  be  to  compare 
the  new  codes  predictions  with  the  results  of  the  physical 
experiments  and  then  to  study  various  free-surface  phenomena. 

The  above  work  is  being  written  up  for  publication  in  a 
series  of  papers.  Tentative  titles  for  these  are  as  follows: 

W.  0.  Pritchard,  Y.  Renardy  and  L.  R.  Scott,  “Tests  of  a 

numerical  method  for  viscous,  incompressible  flows.  I: 
Fixed-domain  problems,"  (in  preparation). 

W.  0.  Pritchard,  Y.  Renardy  and  L.  R.  Scott,  "Tests  of  a 

numerical  method  for  viscous,  incompressible  flows.  II: 
Free-surface  problems,"  (in  preparation ) . 

A  third  publication  is  planned  jointly  by  Malkus,  Pritchard,  Scott 
and  Strikwerda,  but  a  title  has  not  yet  been  chosen. 


Although  the  codes  developed  by  Pritchard  and  Scott  will 
solve  a  wide  range  of  problems,  it  was  also  desired  to  develop 
entirely  new  codes  implementing  new  numerical  techniques.  Two 
independent  directions  were  pursued  during  this  period.  One 
involved  higher-order  spatial  accuracy  and  better # representation 
of  the  incompressibility  constraint  in  the  numerical 
discretisation  based  on  implementing  a  quart ic  finite  element 
method  that  has  been  studied  theoretically  by  Scott  and  Michael 
Vogelius  (of  the  University  of  Maryland).  The  other  addressed 
temporal  discretisation,  especially  with  regard  to 
high-Reynolds-number  flows,  based  on  a  streamline-diffusion 
f inite-element  technique  developed  by  Claes  Johnson  of  the 
University  of  Oothenberg  when  he  was  visiting  the  University  of 


Michigan.  To  achieve  these  goals  most  cost-effectively,  it  was 
decided  to  hire  a  professional  programmer  rather  than  rely  solely 
on  Pritchard  and  Scott  (and  other  collaborators)  for  the  coding,  / 
as  had  been  done  in  the  past.  This  approach  is  quite  standard  in  \ 
governmental  and  industrial  laboratories,  but  it  is  somewhat 
unusual  in  academic  circles.  Although  this  approach  was 
complicated,  primarily  due  to  its  novelty  in  a  university  setting, 
it  was  successful  and  would  be  recommended  for  future  projects 
under  appropriate  circumstances.  (The  availability  of  an 
appropriate  person  is  critical;  our  programmer  left  the  university 
at  the  end  of  this  project  to  work  at  an  industrial  laboratory  at 
a  50k  increase  in  salary  and  would  not  be  enticed  back  to  the 
academic  setting  easily.  Moreover,  the  demand  for  such  people  in 
industry  far  exceeds  the  supply. ) 


The  numerical  studies  with  the  quartic  finite-element  method 
for  viscous,  incompres* * ble ,  two-dimensionsl  flows  were  based  on  s 
modification  of  the  IMSL  code  TWODfcPEP.  Tw^  spprearh**  w»r» 
followed.  The  first,  which  involved  a  straightforward  use  of  the 
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eo 4m,  taitad  the  standard  penalty  method  for  imposing  tha 
incompressibility  eoastraiat.  Typical  results  of  this  study  aro 
fives  in  Tables  1  sad  2.  A  preliminary  description  of  thea  was 
fives  by  8cott  in  tha  conference  paper  "Spatial  discretisation 
techniques  for  the  Mavier-Stokas  aquations:  thaoretieal  and 
coeputational  results, "  which  appeared  in  tha  proceedings  ant it lad 
OOHPOTI NO  METHODS  IN  APPLIED  SCIENCES  AND  ENQINEERING  VI  (R. 
Glowinski  and  J.-L.  Lions,  ads.,  Ansterdan:  North-Holi  and,  .1984 , 
481-470).  Note  that  the  results  in  Table  1  indicate  the  high 
order  of  accuracy  of  the  nethod,  at  least  in  tha  LI  no  re.  (Tha 
poblae  being  solved,  as  shown  in  Figures  9a  and  9b,  is  not 
sufficiently  regular  to  prove  optimal-order  aoouracy,  dua  to  tha 
lack  of  seoothness  of  the  boundary. )  The  data  in  Table  1 
indicates  only  4th  order  convergence  in  the  L2  norm,  and  not  the 
bast-possible  5th  order  convergence  that  approximation  theory 
would  allow.  However,  the  LI  convergence  order  is  about  4.5  at 
tha  last  mash-halving.  .  Also  note  that  the  pressure  accuracy  is 
disappointing,  as  would  be  expected  from  the  penalty  aathod. 

Table  2  indicates  that  reasonable  accuracy  is  retained  even  for 
moderately  high  Reynolds  numbers  and  a  small  number  of  variables. 

Tha  second  approach  involved  a  modified  penalty  method, 
which  holds  the  promise  of  yielding  exact  satisfaction  of  the 
incompressibility  constraint  as  well  as  more  accurate  pressure 
approximations .  Implementing  the  modified  penalty  method, 
together  with  Newton's  method  for  treating  the  nonlinear  term, 
required  a  substantial  reworking  of  the  TffODEFEP  code. 

Convergence  of  the  latter  approach  had  been  cast  in  doubt  in  a 
remark  in  the  book  AUGMENTED  LAGRANQIAN  METHODS  by  Fortin  and 
Glowinski.  and  our  investigations,  as  indicated  in  Table  1,  tended 
to  support  this  observation.  Thus  further  investigation  is  needed 
to  determine  a  method  to  link  the  modified  penalty  approach  with 
efficient  techniques  for  resolving  the  nonlinear  terms. 

The  implementetion  of  the  streamline-diffusion,  finite- 
element  method  for  inviscid  flows  required  starting  a  code  from 
scratch.  The  first  step  was  to  develop  a  solver  for  a  linear, 
hyperbolic  equation  using  the  streamline-diffusion,  finite- 
element  method.  This  is  later  to  be  coupled  with  an  elliptic 
solver  (e.g.,  Randy  Bank's  multi-grid  code)  to  implement  a  full 
nonlinear  Euler  equation  solver.  However,  the  critical  new 
component  is  the  streamline-diffusion  technique  for  temporal 
integration  of  an  advection  equation,  so  this  was  addressed  first. 
Numerical  teats  were  performed  on  this  technique  with  the  code. 

It  appears  that  the  code  is  very  effective  in  limiting  cross-wind 
diffusion  for  discontinuous  solutions  (see  Figure  4),  and  thus 
shows  promise  as  a  component  in  a  high-Reynolds-number  code. 

Teats  were  also  performed  on  problems  with  smooth  data,  cf. 

Figure  5.  Convergence  test  for  this  problem  indicate  seoond -order 
convergence,  and  the  L-2  error  for  the  time  t  s  0.6  shown  in 
Figure  8  was  0.012  .  What  remains  to  be  done  is  the  linking  of 
the  elliptic  solver  as  described  previously  together  with 
development  of  efficient  techniques  for  solving  some  implicit 
equations  that  are  inherent  to  the  nethod  at  each  time  step. 
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TABLI  1 


Converg ence  testa  of  the  quartic  f  ini  to  ol— nt  method  on 
Jaffery-Bamel  flow  os  shown  in  Figure  3.  All  runs  with  ths 
■mm  Reynolds  number:  Rave  =  6.85,  Rux  =  9.88. 


6  tri¬ 

velocity 

velocity 

velocity 

pressure 

»* 

See 

angles 

LI  error 

L2  error 

Max  error 

Max  error 

Kay 

used 

(relative) 

(relative) 

(relative) 

( relative ) 

below 

7 

0.762X-2 

0.131X-1 

0.273E-1 

n/e 

A 

7 

0.67  IE- 2 

0.143E-1 

0.288E-1 

n/a 

B 

7 

0.8711-2 

0.143X-1 

0.288X-1 

0.200X0 

C 

28 

0.6761-3 

0.156X-2 

0.453X-2 

0.108X0 

A 

28 

0. 769E-3 

0.177E-2 

0.508E-2 

0. 122X0 

C 

112 

0.342E-4 

0.117E-3 

0. 566X-3 

0.720X-1 

A 

Kay  to  Ms t hods  Dssd 

A  straight  psnalty,  spsilon  =  lJt-6,  3  stags  continuation,  3 
extra  itsrations  with  full  nonl insari ty 

B  straight  penalty,  spailon  =  1X-6,  3  stags  continuation, 
itsratsd  with  full  nonlinsaritjr  until  MconvsrgoncsH 

C  prsssurs  only  relaxed  (by  1.2***b  at  ths  n-th  itsration) 
psnalty  itsration ,  spsilon  *  11-6,  3  stags  continuation, 
itsratsd  with  full  nonlinsarity  until  "convergence" 


TABLE  2 

Errors  for  ths  quartio  finite  slsnsnt  nethod  calculating 
Jsffsry-HaaMl  flow  for  various  Reynolds  numbers,  using  only 
seven  triangles,  as  shown  in  Figure  3a. 


Reynolds  number  (Rave)  6.9  18  65  222 

Relative  L2  velocity  error  0.13X-1  0.27X-1  0.69E-1  0.15E0 


mesh:  lxl  hole  pressure  17may85 


depths  1,3,2*  1.5  KMESH  l(x2,3extrap)  24sep85 


3ui9^j9av  pug  —  saouaJsyiQ  ssaj^s 


depths  1..75  &  .5  KMESH  I(x2,3extrap)24sep85 


Reynolds  numb 
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